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In nnmerical simnlations stndying preheating in the classical approxima¬ 
tion there is the problem how to derive the classical initial conditions from 
the qnantnm vacnnm flnctnations. In past treatments, the initial conditions 
often pnt an energy density into the classical held of order of the cntoff, lead¬ 
ing to a divergent temperatnre after thermalization. We snggest a solntion to 
the problem which follows natnrally from a Hartree ensemble approximation, 
introdnced recently as an improvement over the standard Hartree approxi¬ 
mation. We stndy the effects on particle nnmbers of the varions treatments, 
within the context of ‘tachyonic preheating’ in 1-1-1 dimensional theory. 
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1 Introduction 


Numerical simulations of quantum field dynamics in real time are currently carried out by 
first making approximations. For systems out of equilibrium the classical approximation 
can be quite useful, for example in studying reheating after inflation In addition 

there are Hartree (gaussian), Dyson-Schwinger and 2PI-functional approximations, possi¬ 
bly supported by large N considerations. These have the advantage of being formulated 
in the quantum theory where potential divergencies are taken care of by renormalization. 

Even if the system is homogeneous it may be important to allow for inhomogeneous 
realizations in order to capture non-perturbative effects. Recently we have formulated a 
Hartree ensemble approximation |P, , in which an overall homogeneous density matrix 

is represented as an ensemble of typically inhomogeneous initial conhgurations, which are 
then evolved in time using the Hartree approximation. With this method we found that 
the poor thermalization properties of the homogeneous Hartree approximation are much 
improved. 

Here we shall describe an application of this method to the problem of spinodal or 
tachyonic instability. This process has been put forward as an efficient mechanism for 
preheating. In ref. i, |] it was studied numerically, using the classical approximation 
with initial values based on quantum fluctuations. However, these initial conditions (of 
the type also used in the other studies in |I|-||^) lead to problems: the hnal temperature 
is of order of the cutoff. Below we propose improved initial conditions for the classical 
approximation which do not have this problem and we compare with results obtained 
using the previously employed initial conditions, in classical and Hartree dynamics. 

We also draw the attention to a fruitful dehnition of particle numbers and energies 
for systems out of equilibrium, which was introduced in |^, and for the fermionic case in 

0 . 

2 Hartree ensemble approximation 

Our tests are carried out in in H-l dimensional theory. The quantum hamiltonian is 
given by 

H = J dx , (1) 

from which follow the Heisenberg operator equations 

(p = TT, TT = (V^ — — A(p^. (2) 

We discretize the held on a space-time lattice, with lattice spacings a and at, such that 
\/(p —>■ {ipx+a — ^x)/ci, 0 —*■ [^xit + at) — (px{t)]/at. The number of spatial lattice sites is iV, 
the volume L = Na and we use periodic boundary conditions. The Hartree approximation 
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implies that the helds can be expanded in terms of mode fnnctions 

^x = ^x + J2 (/“ + fx* bV) , K = Trx + J2 {fx ba + It bi) , (3) 

a a 

in which ip^ and tTx = Px are the mean helds. The creation and annihilation operators 
ha are time-mdependent, with = 5a/3, while the mode fnnctions /" are time- 

dependent. Fnrthermore, the density operator (which is time-independent in the Heisen¬ 
berg pictnre) is ganssian, snch that it is completely specihed by giving the initial valnes of 
the one and two-point fnnctions, {(px) = ^xi {"^x) = t^x, and (pxPy), (px'^y + '^y^x), i'^x^y}- 
With a snitable choice of initial mode fnnctions these can be rephrased as = n^Sap, 

with (ba) = (hi) = (babp) = = 0. In case the ganssian state is pnre, the = 0, 

and the Hartree eqnations for the mean held and mode fnnctions take the form of a set 
of conpled non-linear partial diherential eqnations, 

Ifix = vy, - ip + \ifil + 3A ^ <fii, 

fS = + 3A /;/r) /!'■ (4) 

In 1-1-1 dimensions the mode snm is only logarithmically divergent as the lattice spacing 
a ^ 0, which is taken care of by mass renormalization. In the following we assnme the 
divergence to be snbtracted from the mode snm, snch that /i^ is a renormalized mass 
parameter, with particle mass = p? (symmetric phase) and = —2/r^ (broken 
phase). 

In the Hartree ensemble approximation the density operator p is written as a snper- 
position in terms of ganssian coherent pnre states |i), 

P = (5) 

i 

and the Hartree approximation is applied to each state |z) individnally. Note that p does 
not have to be ganssian or pnre-state. Given p and a choice of the set of coherent states 
|z), the probabilities Pi are nniqnely dehned in simple cases, e.g. a free thermal ensemble 
0]. However, we shall nse (j^) to formulate approximations to the true p, such that the 
initial values of the one and two-point functions are reproduced. We see this as a sort of 
coarse graining: not all details of the true p are kept, only those referring to the important 
low momentum modes of the one and two point functions {px), ■ ■ ■, {^x^y)- 

Even if px = {Px) = J2iPi{i\Px\i) is homogeneous (independent of x), the individual 
(phi = {i\px\i) are typically inhomogeneous. This is important for thermalization, because 
the modes interact non-linearly with the inhomogeneous (phi; the particles can scatter. 
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3 Particle number 


One needs a definition of particle number and frequency ujk which is reasonably intuitive 
and robust, such that it can be applied to systems out of equilibrium. Some coarse graining 
will have to be involved and the following definition takes this to the extreme of using an 
average over the whole volume: we define and Uk at each time by the equations 


(^A:^fc)conn (^k 2^ ^ conn (^k 2^ ^ki 


( 6 ) 


where ^ VT and ^ {‘Pk‘pl) - etc. For a free system 

these formulas produce the standard distribution functions and frequencies, while for 
systems out of equilibrium the above definitions have been used successfully in 1,0- 
A corresponding definition for fermions was introduced in |T^, which also mentions the 


relation to Wigner functions. Here we note in particular that the above equations can be 
easily solved, with positive cok, and we always have found Uk to be positive. 

The above formulas can also be applied to classical fields, with 


5 '^k'^k — 'k^k^k’) 

^k 


(7) 


where ■ ■ ■ denotes a classical average (typically over initial configurations). 


4 Initial conditions 


Consider a quench, a sudden change of the sign of the renormalized mass parameter 
from positive to negative values at time f = 0. This is supposed to mimic the effect 
of a rapid phase transition (cf. figure |l|). Because the potential at t = O’*" is unstable, 
Fourier modes of the field with wave vector \k\ < |/i| are unstable and grow exponentially. 
Occupation numbers grow large and non-linearities become important, which fact has 
been dealt with numerically using the classical approximation. 

The question is how to start the simulation, and, can we do better and stay within 
the quantum theory? Classically, the ground state at t = 0“ implies that the fields and 
canonical momenta are zero, so nothing would happen subsequently. Quantum fluctua¬ 
tions are invoked to start the process. At t = 0“ the system is in a single well ground 
state. Neglecting non-linearities the quantum fluctations are characterized by 


Uk = 0 , ujk = ujI 




( 8 ) 


These expectation values can be reproduced with the classical ensemble distribution 


p(<p,7r) oc n exp 

k 


TTfcTTfc + Ujliflcpk 

(Jk{nk + 1/2) J ’ 


(9) 
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Figure 1: The potential and sketch of the gaussian wave function through a quench. From 
left to right: t = 0“, > 0; t = O'*", < 0; a later time. 


with Ufc = 0 and cvk = The type of initial conditions used in [ll|-|P] correspond 
to generating conhgurations from this distribution and using these as initial values for 
the subsequent classical evolution. A problem with this is that the total energy density 
diverges oc for lattice spacing a —>• 0, and after equilibration the classical 
temperature T also diverges oc a~^. Renormalization of this divergence is not possible out 
of equilibrium (in equilibrium one can attempt to do this with temperature dependent 
countertems BUI)- 


The Hartree approximation does not have this problem as it is formulated in the quan¬ 
tum theory. However, with homogeneous initial conditions it has unsufficient scattering 
to obtain thermalization. 

To improve this situation we proceed as follows. We solve the Heisenberg equations 
analytically for small times, when the modes grow exponentially oc — |/cP). 

Non-linearities become important roughly when = Xip‘^/4, at the spinodal time 

tsp = (l/2|/i|) ln(2|/ip/A). At some resampling time 0 < Rg < tgp we provide a gaus¬ 
sian distribution (|]), for initial values (pk and tta, of the unstable modes, based on the 
analytically calculated two-point functions {(fl^Pk) and {^\'^k) (in the analytic solution 
{^k) = (TTfc) = {^l^k + ^k^k) = 0): 


- In the classical approximation, the distribution is chosen of the form (|^) with 
Ufc -|- 1/2 —> Uk (so dropping the 1/2), such that the analytically calculated par¬ 
ticle distribution Uk is reproduced by (|^). The high momentum modes in (^ are 
suppressed by nk dropping to zero (like as it turns out). We actually keep only 
the unstable modes in the initialization as for these typically S> 1. The resulting 
energy density is hnite. 


- In the Hartree ensemble approximation we also use the form (P) to generate initial 
values for the mean helds in each Hartree realization, again with -|- 1/2 —> n^. 
The mode functions are initialized as vacuum plane waves a ^ k, 

fx = / + F, = -2/i2, (10) 

pLuT 

where m is the mass of the particles and the gaussian coherent states are specihed 
by hk\i) = 0. (It would perhaps have been logically more consistent to use uik in 
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Log( 1+1/n) tm = 10 Log( 1+1/n) tm = 2980 




Figure 2: Particle distributions at various times using ‘old’ initial conditions. The straight 
line is a £t to the data at small ujk- The classical distribution rik = T/ujk with the same 
temperature is also indicated (curved line). 

stead of in which case the quantum two point functions at t = trs would have 
been reproduced exactly, with the l/2ci;fc in provided by the contribution of the 
mode functions.) Since the mean helds in each realization are inhomogeneous we 
expect much improved thermalization. 


5 Numerical tests 

We used relatively weak coupling, = rn?/2\ = 6, a lattice spacing am = 0.2 and volume 
Lm = 102.4. Figure ^ shows the particle distribution for several times, obtained by 
classical evolution using the ‘old’ initial conditions at t = 0+, i.e. (P) with nfe + 1/2 —>■ 1/2 
(so keeping the 1/2) and ojk = a;°. Plotted is ln(l + l/n/) versus which would give 
a straight line with slope 1/T for a Bose-Einstein distribution = exp(a;fe/T) — 1. 
One clearly sees at time tm = 10 that the high momentum modes, which have not been 
affected yet by the evolution, are populated according to the vacuum fluctuations: their 
effective classical particle number (from (|]) and (|^) is 1/2 and ln(l + 1/0.5) ~ 1.1. At 
later times the distribution approaches the classical T/co^. 

Next we show results using the proposed ‘new’ initial conditions for the classical evo¬ 
lution. With the parameters used here there are then 23 unstable modes (12 with k >0). 
The spinodal time t^^m ~ 1.8 and the resampling time was chosen to be t^-gm ~ 1. Figure 
H shows the results for this case. We see that initially Uk drops very fast as a function of 
k (cf. tm = 10). At intermediate times the distribution resembles a Bose-Einstein form. 
The approach to the classical distribution is much slower than in the previous case, the 
flow of energy to the higher momentum modes is evidently slow. 

Results using the Hartree ensemble approximation are shown in hgure ^left. In this 
case we see a slow approach to a Bose-Einstein distribution, within the same time span 
as in the previous cases. However, we expect the distribution to approach a classical 
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Log( 1+1/ n) tm = 10 



Log( 1+1/ n) tm = 2830 



Figure 3: As in figure using the new initial conditions. 



Figure 4: Left: as in figure ^ using the Hartree ensemble approximation. Right: temper¬ 
atures as a function of time, for the classical evolution with ‘old’ initial conditions (top), 
the classical evolution with new initial conditions (bottom) and the Hartree ensemble 
approximation (middle). 


form on a much larger time scale, typically by a factor of 100 . 

Finally, figure fright shows how the temperature (as determined by the slope in figs, 
at low momenta) evolves in time for the three cases. In the first case we see the 
temperature gradually rising, eventually T will get of order 1/a. In the second case T 
drops slowly as the energy is being equipartitioned from the low momentum modes to all 
modes. For the Hartree ensemble case the temperature remains reasonably constant after 
the initial transient. We checked that the same temperature is obtained using a smaller 
resampling time. 


6 Conclusion 

The divergences in the energy density and the classical temperature are cured by our pro¬ 
posal for obtaining initial conditions for classical simulations. It makes a big difference at 


6 



















later times. The Hartree ensemble approximation appears to be a significant improvement 
over the classical treatments. However, this method is nnmerically expensive and fnrther 
tests are needed in 3+1 dimensions. 
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